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1. Introduction - e + e Colliders 

The production of hadronic final-states by a variety of interactions is illustrated 
in Fig. 1. In electron-positron annihilation hadronic activity is, by construc- 
tion, limited to the final state, making the study of hadronic events cleaner and 
simpler relative to lepton-hadron and hadron-hadron collisions, from both the ex- 
perimental and theoretical points-of-view. On the experimental side there are no 
remnants of the beam particles to add confusion to the interpretation of hadronic 
structures, and, apart from initial and final-state photon radiation effects, the 
hadronic centre-of-mass frame coincides with the laboratory frame. On the the- 
oretical side the absence of hadrons in the incoming beams removes dependence 
on the limited knowledge of the parton density functions of hadrons, as well as 
rendering QCD calculations at a given order of perturbation theory easier to per- 
form because there are generally fewer strong-interaction Feynman diagrams to 
consider. Electron-positron annihilation thus provides an ideal environment for 
precise tests of QCD. 

A large number of e + e~ colliders have been constructed over the past 25 years; 
these are listed in Table 1. The range of cm. energies Q extends from a few 
GeV at the very first colliders up to almost 200 GeV at the CERN LEP-II col- 
lider. The first generation of colliders was built on speculation of allowing exciting 
high-energy physics studies. They did not disappoint, the J/\& being discovered 
at SPEAR, the gluon being observed at PETRA, and a wealth of strong- and 
electroweak-interaction studies being performed at PETRA, PEP, and TRISTAN, 
all of which served to establish the validity of the Standard Model. A second gen- 
eration of colliders has been designed to serve as particle 'factories': DA$NE in 
the vicinity of the $ resonance, BEPC near the charmonium threshold, DORIS, 
CESR, PEP-II and KEKB around the bottomonium resonances, SLC and LEP 
at the Z° resonance, and LEP-II at the W + W" threshold. 
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Figure 1: Schematic of the production of hadronic final-states by different interac- 
tions: e + e~ annihilation, deep inelastic lepton-hadron scattering, hadron-hadron 
collisions, and production and decay of heavy quarkonia. 
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Collider 


Location 


cm. energy Q (GeV) 


ADONE 


Frascati 


1 - 3 


DCI 


Orsay 


1 - 2.4 


• VEPP-2M 


Novosibirsk 


1 - 1.5 


• DA$NE 


Frascati 


1 - 1.5 


SPEAR 


SLAC 


2-8 


• BEPC 


Beijing 


2-3 


DORIS 


DESY 


3-11 


• VEPP-4M 


Novosibirsk 


10 - 12 


• CESR 


Cornell 


10 - 11 


• PEP-II 


SLAC 


9(g) 3.1 


• KEKB 


KEK 


8 3.5 


PETRA 


DESY 


12 - 47 


PEP 


SLAC 


29 


TRISTAN 


KEK 


50 - 64 


• SLC 


SLAC 


88 - 93 


LEP 


CERN 


88 - 93 


• LEP-II 


CERN 


130 - 192 


XLC 


???? 


500 - 1500 



Table 1: e + e colliders 1972 - 200? • denotes running/under construction. 

With the exception of the SLC, all of these colliders have been of the storage 
ring type, the largest, LEP-II with a circumference of 27km, probably marking the 
limit of the energy that can be achieved with current storage ring technology for an 
acceptable cost. The SLC is the first example of a high-energy linear e + e~ collider; 
it achieves the same collision energy as LEP, but has an effective length of about 
3 miles and was considerably cheaper to construct. Because of their intrinsically 
lower cost/GeV, linear colliders represent the obvious path towards construction 
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of higher-energy e + e~ colliders with current acceleration technology. A number 
of proposals for such an accelerator are represented by 'XLC in Table 1; they 
all aim to achieve cm. energies between 500 and 1500 GeV, which is believed to 
cover the interesting range for study of electroweak-symmetry-breaking processes. 
Some examples of QCD tests that could be made at the XLC will be given towards 
the end of these lectures. 

It would require a semester-long lecture series to do full justice to QCD studies 
in e + e~ annihilation, so some hard choices have been made as to the material to be 
covered here; I apologise well in advance for all that has been omitted. No attempt 
has been made to give a complete review of all of the experimental results in any 
of the areas covered; usually one or two results or figures are shown as examples. 
For this purpose I have drawn heavily on material from TASSO and SLD, the two 
experiments with which I have been involved since 1985; no disrespect is intended 
to the many other experiments whose results may not be shown. Tests of QCD 
in hadron-hadron and lepton-hadron collisions will not be discussed here as they 
are covered in other lectures jl], |2| at this Institute. 

In the interests of pedagogy I shall review the fundamental properties of QCD 
and the important experimental measurements from e + e~ annihilation that have 
been key historically to establishing the theory. Having verified QCD, in a qual- 
itative sense, as being the only viable theory for describing strong interactions, I 
shall then review quantitative tests in the form of measurements of a s , the single 
parameter of the theory, and put the e + e~ measurements into context with deter- 
minations from other processes. I shall focus on measures of the event topology, 
especially on jet definition and the relation between the jets observed in detectors 
and the 'true' underlying parton-jet structure. This will introduce the problem of 
hadronisation, as well as the difficulty of relating finite-order perturbative QCD 
calculations to the 'all-orders' data. I shall conclude by looking forward to the 
precise QCD tests that could be made at a high-energy e + e~ collider and to the 
qualitatively new ttg system accessible at such a facility. 
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2. Our Theory of Strong Interactions - QCD 



Quantum Chromodynamics (QCD) || is our theory of the strong interaction 
between quarks and gluons. It is a non-Abelian Yang-Mills gauge theory that 
describes the interactions of a triplet of spin-1/2 quarks possessing the colour 
quantum number (c = r,b,g) via an octet of vector gluons. The spinor quark fields 
q c {x) transform as the fundamental representation of the SU(3) group, whilst the 
gluon fields A a ^(x) (a = 1,2,. . .,8) transform according to the adjoint representa- 
tion. The SU(3) colour transformations are generated by the 3x3 matrices T a 
= X a /2, where A a are the Gell-Mann matrices |f] which obey the commutation 
relations: 

[T a ,T b ] = if abc T c (1) 

and f abc are the structure constants of SU(3). The Lagrangian has the form: 

C = -\F; v F^ a + q (ij^ - m) q (2) 

where F^ u is the field strength tensor: 

F% = d»Al-d v Al + gr hc A\Al (3) 

and _D M is the covariant derivative: 

D, = d,-igT a Al{x), (4) 

g is the bare coupling of the theory, m the bare mass of the quark field and the 
gluons are massless. 

Following ||, the 'essential features' of QCD may be summarised as: 

• quarks with spin 1/2 exist as colour triplets 

• gluons with spin 1 exist as colour octets 

• the coupling qqg exists 

• the couplings ggg and gggg exist 
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• the couplings are equal 

• the coupling decreases as 1/lnQ 2 

For most of the first lecture I shall review the evidence, from e + e~ annihilation alone, 
that QCD is the correct theory of strong interactions. 

3. Establishing the QCD Lagrangian 
3.1 Two- Jet Events and qq Production 

Quarks were first postulated in 1964 by Gell-Mann and Zweig || as a calculational 
device to explain the rich spectroscopy of recently-discovered mesons and baryons 
in terms of bound qq and qqq (or qqq) states. The first direct evidence for quarks 
came from the observations at SLAC in the late 1960s that in electron-nucleon 
scattering experiments at high Q 2 the electron scatters from quasi-free pointlike 
particles. In e + e~ annihilation a convincing demonstration of the existence of 
quarks was provided by the observation of jets in the Mark I experiment at SPEAR 
in 1975 This analysis represents the first use of an event shape observable 
which, as will be discussed later, is a key component in the study of hadronic final 
states, so I shall briefly describe it. 

By eye the spatial distribution of particles in hadronic events recorded in the 
Mark I detector operating at cm. energies between 3.0 and 7.4 GeV looked 
more-or-less isotropic, and it was hard to distinguish any clear jet structure. The 
quantity sphericity, 

S = M " (S t i - ) , (5) 

where Pi represents the momentum of particle i and the sums run over all particles 
in each event, was invented || to characterise the degree of isotropy in the particle 
flow. In each event an axis, the sphericity axis, is defined so as to minimise the 
quantity in brackets in the numerator; eq. (5) then defines the sphericity of the 
event. A completely isotropic distribution of particles, or spherical event, would 
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yield S ~ 1, whilst a perfectly-collimated back-to-back two-jet event would have 
S — 0. Sphericity distributions from Mark I are shown in Fig. 2 for data taken at 
several different cm. energies. As the energy was raised from 3.0 to 7.4 GeV a 
clear change in the sphericity distribution was observed, the distribution shifting 
to lower values at higher energies. This was interpreted in terms of an increasing 
degree of collimation of particle production with cm. energy, namely the onset 
of the production of two back-to-back jets of hadrons. At higher energies the jet 
structure is much more apparent by eye, as indicated in the Z° decay event from 
SLD shown in Fig. 3, and is striking evidence for the production of a back-to-back 
quark and antiquark in e + e~ annihilation. 

The Mark I analysis was also able to establish the nature of the spin of the 
quark and antiquark. Shown in Fig. 4 is the distribution of the azimuthal-angle, 
<f), of the sphericity axis w.r.t. the beamline, at two cm. energies. At Q = 7.4 
GeV the electron and positron beams in the SPEAR ring built up a degree of 
transverse polarisation P via the Sokolov-Ternov synchrotron radiation effect || 
and a clear modulation in <fi is visible. This is in contrast to the flat </> distribution 
at Q = 6.2 GeV which corresponds to a beam-depolarising resonance (P = 0) in 
the accelerator. A fit of the function: 

— oc 1 + ol cos 9 + P asm 9 cos 20 (6) 

to the 7.4 GeV data yielded a = 0.78±0.12; this is close to unity, which is expected 
for production of two spin-1/2 particles flOf . 



S 



80 




Figure 2: Sphericity distributions measured by the Mark I experiment at 
SPEAR at cm. energies of (a) 3.0, (b) 6.2, (c) 7.4 GeV. The narrowing 
of the distribution, and the trend towards smaller values as the cm. energy is 
raised, represent evidence for collimated production of hadrons in e + e~ annihila- 
tion. The dashed line represents the expectation from a 'phase-space model' of 
hadron production, (d) As (c) but for a subset of events containing particles with 
scaled momentum, 2p/Q, less than 0.4. 
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Run 10507, EVENT 346 

27-MAR-1992 06:59 

Source: Run Data Pol: 

Trigger: Energy CDC 

Beam Crossing 1119271 




Figure 3: A contemporary two-jet event recorded by SLD: the process 
e + e~ — > Z° — > qq. 

These studies were subsequently extended at the higher-energy PETRA col- 
lider, and examples from TASSO |l l[ are shown in Fig. 5. Here the distribution 
of the polar-angle (9s) of the sphericity axis is shown at cm. energies of 14, 22 
and 35 GeV. A fit to the functional form: 

f N a oc 1 + a S)T cos 2 6 s ,t (7) 
a cos 

yields, at 35 GeV for example, as = 1.03±0.07, again characteristic of the produc- 
tion of two spin-1/2 particles in the e + e~ annihilation. Also shown in Fig. 5 is our 
second example of an event shape observable in the form of the thrust-axis |p~2f 
polar-angle (6*t) distribution. Thrust will be discussed later; it is qualitatively 
similar to sphericity in that it can be used to quantify the degree of collimation of 
particle production, although it has properties that make it more attractive the- 
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Figure 4: Azimuthal-angle distribution of the sphericity axis from Mark 10. At 
Q = 6.2 GeV (a) the SPEAR beams were unpolarised and at Q = 7.4 GeV (b) 
the polarisation-product was 0.47; the modulation in (b) is characteristic of the 
production of two spin-1/2 particles in e + e~ annihilation. 

oretically. The thrust-axis polar-angle distribution in Fig. 5 was fitted to obtain, 
at 35 GeV for example, a T = 1.01 ±0.06, in good agreement with the result using 
the sphericity axis. 

So far spin-1/2 quarks and antiquarks would appear to be well established, 
and their colour-triplet nature, Nq = 3, is required in the quark-parton model 
(QPM) of hadrons to explain the existence of spin-3/2 baryon states such as the 
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Figure 5: Polar- angle distributions of the sphericity and thrust axes from 
TASSO [|ll]]. The 1 + cos 2 9 form is characteristic of the production of two spin-1/2 
particles in the e + e~ annihilation. 

A ++ (u|u|u|) and Q~ (s|s|sf), which would otherwise contain three identical 
fermions in the same quantum state, in violation of the Pauli exclusion principle. 
In e + e~ annihilation evidence for Nq = 3 is provided by the quantity: 

_ cr(e + e~ — > hadrons) 
erQ£Z)(e + e- _> /j+fj,-) 

which, according to QED and the QPM, should be equal to iVc E/ qj, where qf 
is the charge of the quark of flavour / and the sum runs over all active flavours 
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at a given cm. energy. 
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Figure 6: The R ratio as a function of cm. energy ||13|| . The expectation for Nc 
= 3 is shown as the solid line. 

A summary of R measurements made up to 1988, as a function of cm. en- 



ergy, is presented in Fig. 6 [0. This is a tremendously information-rich figure. 
First, increases in R just above Q 2 = 10 and 100 GeV 2 represent the cc and 
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bb production thresholds - further evidence, were it needed, for the existence of 
quarks. Secondly, the QED + QPM prediction comes close to the data only if 
the quarks are assigned fractional charges and the number of colours Nc = 3 is 
used; the 'colour singlet' expectation (Nc = 1) is simply too low by a factor of 
about three! Thirdly, above Q 2 = 1000 GeV 2 the data points rise as Q 2 increases, 
representing the onset of contributions to e + e~ annihilation from Z° exchange. 
Finally, in regions between quark flavour thresholds and below the tail of the Z° 
resonance, there is a residual excess in the data relative to the QED + QPM 
expectation, and the excess appears to decrease as Q 2 increases. In other words, 
some mechanism causes an increase in the 'phase-space' for hadron production 
beyond QED + QPM, but at a rate that decreases with Q 2 . In the language of 
the 1990s we know that the extra contribution is due to gluon emission in the 
final state, and that the probability for this process, a s , decreases roughly loga- 
rithmically with Q 2 . The i?-ratio thus provides indirect evidence for the existence 
of the gluon, as well as for the non-Abelian 'running' of the strong coupling. 



3.2 Three-Jet Events and the Gluon 

In e + e~ annihilation events containing three distinct jets of hadrons were first 



observed in 1979 at the PETRA storage ring at cm. energies around 20 GeV. 
Such events were interpreted [|15| in terms of the fundamental process e + e~ — > qqg, 



providing direct evidence for the existence of the gluon and its coupling to quarks. 
A modern example of a three-jet event, in fact the very event used to advertise 
this Summer Institute, is shown in Fig. 7. 

Counting the number of jets per event, and then comparing the numbers of 



two- and three-jet events, it was found |T6[ that around Q = 20 GeV 

t 3 ~ jet6VentS - 0.15. (8) 
#2 — jet events 

Since, at lowest order in perturbative QCD, this ratio is simply the probability 
for gluon emission, or a s , the strong coupling parameter, simple event counting 

14 




Figure 7: A contemporary three-jet event recorded by SLD: the process 
e + e~ — > Z° — > qqg. 

indicated that the strong coupling at 20 GeV was around 0.15, i.e. about ten times 
larger than the electromagnetic coupling a. More systematic determinations of 
a s will be discussed later. 

Having observed the gluon directly in three-jet events one still needs to know 
whether it is the gluon of QCD, namely a colour-octet vector particle. Many 
studies of the nature of the gluon spin were performed at the PETRA and PEP 
storage rings and involved analysis of the partition of energy among the three jets. 
Ordering the three jets in e + e~ — > qqg according to their energies E\ > E 2 > E 3 , 
and normalising by the cm. energy Q, we obtain the scaled jet energies 

Xi = ^ (2 = 1,2,3), (9) 

represented in Fig. 8, where x\ + x 2 + £3 = 2. Making a Lorentz boost of the 
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event into the rest frame of jets 2 and 3 the historically-important Ellis-Karliner 



angle 9ek is defined []T7| to be the angle between jets 1 and 2 in this frame. For 
massless partons at tree- level: 

cosQ EK = (10) 

Xi 

The results of an early study by TASSO are shown in Fig. 9, where the Ellis- 
Karliner angle distribution is compared, for data taken at Q ~ 30 GeV, with 
the prediction of QCD. One can also consider alternative 'toy' models of strong 
interactions, for example a model incorporating spin-0 (scalar) gluons [[0|] . From 
Fig. 9 the scalar-gluon model is clearly excluded. 





Figure 8: (a) Representation of the momentum vectors in a three-jet event, and 
(b) definition of the Ellis-Karliner angle. 

Similar studies have been extended at the Z° resonance by the LEP and SLC 
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Figure 9: The Ellis-Karliner angle distribution of three-jet events recorded by 



TASSO at Q ~ 30 GeV JTS); the data favour spin-1 (vector) gluons. 



experiments. In this case the inclusive differential cross sections, calculated at 
leading order and assuming massless partons, can be written for vector gluons 



1 d 2 a v 



oc 



x\ + x\ + (2 — x\ — X2) 



a dx\dx2 (1 — xi)(l — x 2 ){x\ + x 2 — 1) 



for scalar gluons |19| : 



1 d 2 a s 
a dx\dx 2 

where 



xf(l - xi) + - x 2 ) + (2 - x 1 - x 2 ) 2 {xi + x 2 - I) 
(1 - xi)(l - x 2 ){xi + x 2 — i) 



-R 



R 



10 Ej-aJ 
Zj(v] + of) 



(12) 



(13) 
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and ctj and Vj are the axial and vector couplings, respectively, of quark flavor j 
to the Z°, and for a model of strong interactions incorporating spin-2 (tensor) 



gluons [£C], [21 



1 d 2 a T (xi + a;2-l) 3 + (l-xi) 3 + (l 
oc - 



X2 



(14) 



a dxidx 2 (1 — X\){1 — x 2 )(xi + x 2 — 1) 

Singly-differential cross sections for x-y, x 2 , x% or cosOek can be obtained by numer- 



ical integrations of Eqs. (11), (12) and (14) and are compared with SLD data |21 
in Fig. 10. The shapes are different for the vector, scalar and tensor gluon cases 
and only the vector case describes the data. 
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Figure 10: Comparison of leading-order scalar-, vector-, and tensor-gluon models 
with SLD data pi| ; the scalar and tensor hypotheses are clearly excluded. 



An additional interesting observable in three-jet events is the orientation of 
the event plane w.r.t. the beam direction, which can be described by three Euler 
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angles (Fig. 11). These angular distributions were studied first by TASSO [ f2"2"f , 
and more recently by L3 p3j and DELPHI |24fl . Again, the data were compared 



with the predictions of perturbative QCD and a scalar gluon model, but the Euler 
angles are less sensitive than the jet energy distributions to the differences between 



the two cases ||23|| . One can parametrise the angular distributions in the form: 

da 



dcos9 
da 



cfcos^Tv 



oc 1 + a(T)cos 2 6 (15) 
oc 1 + ajv(T)cos 2 #jv (16) 



^ oc 1 + /3(T)cos2 X (17) 
where T is the thrust value |l2j of the event. As an example, fits of eq. (16) 



to SLD distributions of cosOjy are shown in Fig. 12 |^TJ in four bins of thrust. 
The coefficients a(T), «jv(T) and (3{T) depend on the gluon spin; they are shown 
in Fig. 13 for leading-order calculations incorporating vector, scalar and tensor 
gluons. The measured a(T), ajq{T) and (3{T) are also shown in Fig. 13 and 
confirm that only vector gluons are compatible with the data. 

At this point it is worth pausing to take stock of what has been learned so 
far. The e + e~ — > two-jet events have provided direct evidence for qq production, 
and the jet axis angular distribution indicates that the quark and antiquark have 
spin-1/2. From the inclusive .R-ratio we can confirm the fractional nature of the 
quark charges, and know that the quarks must exist as colour triplets since Nc = 
3 is the only value that brings QED + the quark-parton model close to the data. 
The value of the -R-ratio also tells us that there must be contributions to hadronic 
final states in addition to qq production, and we know that these are provided 
by three-jet events, which represent direct evidence for the existence of the gluon 
and its coupling to quarks and antiquarks. The distributions of jet energies, or 
equivalently of jet angles within the event plane, as well as of the event plane 
orientation itself, confirm that the only hypothesis that fits the data is that the 

19 




Figure 11: Definition of event-plane orientation angles. 

gluon has spin-1, and therefore that it is the vector boson of QCD. Finally, from 
counting the relative rates of three- and two-jet events at Q ~ 30 GeV we know 
that the coupling strength of the gluon to quarks is about 0.15. Checking the list 
of 'essential features' of QCD we see that we have verified about half of them! 
The next item in the list refers to the triple- and quartic-gluon couplings; in order 
to study these we need to examine multi-jet final-states. 
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Figure 12: Distributions of cosO^ as a function of event thrust from SLD |plj| 
3.3 Multi-Jet Events and Gluon Self- Couplings 



Consider the Feynman diagrams for production of 4-jet final states shown in 
Fig. 14. Figs. 14(a) and (b) illustrate the gluon Bremsstrahlung process, whilst 
Fig. 14(d) shows the splitting of a gluon into a quark and an antiquark; the latter 
may be thought of as a QCD analogue of the QED process whereby a photon 
converts into an electron and a positron. Figs. 14(a,b,d) are sometimes referred 
to as 'Abelian' diagrams. Fig. 14(c) illustrates the lowest-order diagram for a 
gluon to split into two gluons. This process has no analogue in QED since the 
photon does not couple to itself, and is a consequence of the non- Abelian nature 
of QCD in that the gluons, by virtue of possessing colour charge, can interact 
among themselves. 
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Figure 13: Dependence of the coefficients of the three-jet event plane orientation 



angular distributions on event thrust value [21 



Now consider the formal properties of the SU(3) group. The group can be 
characterised by constants known as Casimir factors that are defined by: 
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Figure 14: Tree-level Feynman diagrams for 4-jet production in e + e annihilation. 

Z a (T a T a ) tJ = 5 tJ C F (18) 

^a, b (f abc f abd ) = S cd N c (19) 

Tr[T a T b ] = 5 ab T F (20) 

The Casimir factors for several common groups are shown in Table 2. We see that 
in the case of SU(3) N c corresponds to the now-familiar 'number of colours' that 
we have already encountered several times. The tree-level couplings appearing 
in Fig. 14 may be classified in terms of the Casimir factors, as illustrated in 
Fig. 15. The amplitude-squared corresponding to the Bremsstrahlung diagrams 
(Fig. 14(a,b)) is proportional to C F , that corresponding to g — > qq is proportional 
to T F , and that corresponding to the non-Abelian process g — * gg is proportional 
to N c . 
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Group 


N c 


c F 


Tf 


U(l) 





1 


1 


U(l)s 





1 


3 


SU(N) 


N 


{N 2 - 1)/2N 


1/2 


SU(3) 


3 


4/3 


1/2 



Table 2: Casimir factors for some common gauge groups. 
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Figure 15: Casimir classification of tree-level QCD couplings. 

It is interesting to consider whether the Casimir factors of SU(3) QCD can be 
measured. Clearly nature does not deliver events corresponding to the tree-level 
vertices shown in Fig. 15! Instead, one must write down the Feynman amplitudes 
for the 4- jet event diagrams shown in Fig. 14, add them to those for 2- and 3-jet 
production at the same order of perturbation theory, and square them to derive 
the total hadronic cross section. The terms corresponding to 4- jet production 
can then be identified in a gauge-invariant manner, and yield a differential cross 
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section of the form: 



1,4 (ol s C f 
— da = 1 



C V 7T 



+ 



a 8 C F s 1 
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2r / IN C \ N c 



^N,F D+ U-\^)F E 



(21) 



where -F4 • • • -^b are kinematical functions. We see that the overall normalisation 
of the cross-section is proportional to (a s Cf) 2 , and that the kinematical distribu- 
tion of the four jets depends on the ratios Nc/Cp and Tf/Cf, which can hence 
in principle be measured. 

The issue of jet definition will be discussed in detail in Section 4.4. For now 
let us assume that 4-jet events can be defined and measured in particle detec- 
tors, and that they can be related meaningfully to the underlying 4-jet parton 
structure described by eq. (21). Two important physical characteristics underly 
the definition of 4-jet observables that are sensitive to Nq/Cf and Tf/Cf'- the 
first is that the two jets resulting from the primary quark and antiquark produced 
in the Z° decay tend to be more energetic than the jets produced by the two 
radiated gluons or the radiated qq; the second is that in the 'non-Abelian' process 
(Fig. 14c) the two gluons tend to be produced in the plane of the primary quark 
and antiquark, whereas in the 'Abelian' process (Fig. 14d) the radiated quark and 
antiquark tend to be produced along an axis normal to this plane 



With this in mind, a number of 4-jet observables that are potentially sensitive 
to the ratios of Casimir factors have been proposed over the years. If one orders 
and labels the four jets in an event in terms of their momenta (or energies) such 
that pi > P2 > P3 > Pi one can define the Bengtsson-Zerwas angle p6[ (Fig. 16): 

cosxbz oc (pi x p 2 ) ■ (p 3 x p A ) (22) 

and the Nachtmann-Reiter angle (Fig. 17): 

cos6* NR oc (pi - p 2 ) ■ (p 3 - pk)- (23) 
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The sensitivity of these observables is illustrated in Fig. 18, where the distributions 
of these angles are shown for SU(3) QCD, as well as for a straw-person U(l)3 
Abelian model of strong interactions, and are compared with L3 data ||27|| ; the 
Abelian model is clearly excluded. 
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Figure 16: Definition of the Bengtsson-Zerwas angle. 



A more recent analysis by OPAL |28| is summarised in Fig. 19; here a simul- 
taneous fit was performed to the Nachtmann-Reiter and Bengtsson-Zerwas angle 
distributions, as well as to the angle a 34 between jets 3 and 4. The resulting 
values of N c /Cp and TpjCp are displayed in Fig. 20, where they are compared 
with the expectations from numerous gauge groups. The SU(3) QCD expectation 
is clearly in good agreement with the data. The expectations from several other 
gauge models, such as SU(4), Sp(4) and Sp(6), also appear to be compatible with 
the experimental results. Note, however, that none of these models contains three 
colour degrees of freedom for quarks, and hence all can be ruled out on that ba- 
sis. Besides SU(3), only the U(l) 3 and SO(3) models contain three quark colours, 
but both are inconsistent with the measured values of Nc/Cp and Tp/Cp. The 
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Figure 17: Definition of the Nachtmann-Reiter angle. 

results shown in Fig. 20 hence yield the remarkable conclusion that SU(3) is the 
only known viable gauge model for strong interactions. 

Recalling the 'essential feature' of QCD that the ggg vertex must exist, we see 
from Fig. 20 that the non-zero measured value of N c /Cf provides direct evidence 
for its contribution to 4-jet production. Now consider the existence of the gggg 
vertex; it should come as no surprise that we need to study events of yet higher 
jet multiplicity in order to be sensitive to it. The tree-level Feynman diagrams for 
5-jet production in e + e~ annihilation are shown in Fig. 21; the gggg vertex can 
be seen in the two diagrams just left of centre on the bottom row. 

Performing a similar exercise to that for the 4-jet cross section one finds: 



L da s = Ida 2 * 3 * + -da 4 ** 1 *, 



(24) 



a (To a 

The first term contributes about 85% of the 5-jet cross section, and may be written: 



-da 



2q3g 



a s C F 

7T 



Gj 



N °-Gs + (^)Gc 



(25) 



where Ga, Gb and Gc are kinematical functions. The contribution of the gggg ver- 
tex is represented by the last term in eq. (25), which is proportional to {Nc/Cf) 2 . 
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Figure 18: Illustration of the sensitivity of the Nachtmann-Reiter and Bengtsson- 



Zerwas angles to non-Abelian effects and comparison with L3 data [27|. 



We have just seen that Nc/Cf must be non- vanishing in order to describe the 4- 
jet data, so that the existence of the gggg vertex is absolutely required in QCD in 
order for the theory to be gauge- invariant and self-consistent. Pushing pedagogy 
to its limits, however, one can still ask if the data actually require the existence 
of the gggg vertex, from a phenomenological point-of-view. One can therefore 
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Figure 19: Simultaneous fit of QCD to OPAL measurements of Xbz, &nr and 
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define a set of ad hoc 5-jet correlation observables, such as those illustrated in 
The measured distributions of the five of these observables that 



Fig. 22 



are most sensitive to the (Nc/Cp) 2 term are shown in Fig. 23, from the OPAL 



Collaboration [29 . 
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Figure 20: The TpjCp vs. N c /Cp plane showing the fitted values derived from 
Fig. 19, as well as the expectations from numerous gauge groups |28[1 ; QCD is in 
good agreement with the data. 



Two possible strategies now present themselves for testing the self-consistencey 
of QCD. One could fit inclusively the quantity Nc/Cf to the 5-jet data shown in 
Fig. 23 and compare it with the value determined from 4-jet events; the results of 
such a comparison are shown in Fig. 24a; the 4-jet and 5-jet events clearly yield 
consistent results. A second possibility is to fit phenomenologically only the gggg 
contribution proportional to (Nc/Cf) 2 '-, the results are shown in Fig. 24b. In the 
latter case the error bars are large due to the small number of 5-jet events, as well 
as to the large uncertainties on multijet production that arise from hadronisation 
effects (see Section 4.4). The measured value of (Nc/Cf) 2 is clearly consistent 
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Figure 21: Tree- level Feynman diagrams for 5-jet production in e + e annihilation. 

with the QCD expectation of (9/4) 2 5, but it is also consistent with zero, 
so that the existence of the gggg vertex has not yet been established from a 
phenomenological point-of-view. 
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Figure 22: Illustration of correlation observables among the jets in 5-jet 
events [p9| . 



3.4 Review of Strategy for QCD Tests 

At this point we have seen that most of the 'essential features' of QCD have been 
established empirically, with the possible exception of the gggg coupling. Even in 
this case, given the existence of the ggg vertex, the gggg vertex must exist in QCD 
in order for the theory to be gauge- invariant. The last 20 years of hadronic-event 
studies at e + e~ colliders have hence established, in a qualitative sense, that the 
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Figure 23: The five correlation observables in 5-jet events that are most sensitive 



to contributions from the gggg vertex [29 



QCD Lagrangian is the correct one to describe strong interactions. At this point 
it therefore seems sensible to revise the strategy for testing QCD. 

Since the theory contains in principle only one free parameter, the strong cou- 
pling a s , QCD can be tested in a quantitative fashion by measuring a s in different 
processes and at different hard scales Q. The precision of these measurements, 
and the resulting degree of consistency among them, determine quantitatively the 
precision with which the theory has been tested. This philosophy is directly anal- 
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Figure 24: Measurements of (a) (N c /Cp) and (b) (N c /Cp) 2 made using 5-jet 



events at OPAL p9| 



ogous to that used to test the electroweak theory by measuring a large number of 
observables that are sensitive to a few key unknown parameters of the theory. In 
addition to testing QCD, the precise measurement of a s allows constraints on pos- 



sible extensions to the Standard Model (SM) of elementary particles; see eg. [[30 



Measurements of a s have been performed in e + e~ annihilation, hadron-hadron col- 
lisions, and deep-inelastic lepton-hadron scattering, covering a range of Q 2 from 
roughly 1 to 10 5 GeV 2 . In the next section I shall describe the e + e~ measurements, 
and compare them with those made in other hard processes; for a review of this 



field see f 3 1 . 
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4. Measurements of a s in e + e Annihilation 

4.1 Theoretical Considerations 

An inclusive observable X may be written schematically: 

X = X EW (1 + 6 QCD ) (26) 

where X EW represents the electroweak contribution. Since, with observables of 
this type, ct s enters via the small QCD radiative correction, 5® CD , a precise mea- 
surement of a s generally requires a large data sample. Observables can also be 
defined that are directly proportional to 5® and hence potentially more sensitive 
to a s . In either case S QCD can be separated into perturbative and non-perturbative 
contributions: 

^QCD _ £pert _|_ ^non-pert (27) 

The perturbative contribution can in principle be calculated as a power series in 
a s , though in practice the large number of Feynman diagrams involved renders a 
complete calculation beyond the first few orders intractable. The non-perturbative 
contribution, often called a 'hadronisation correction' in e + e~ annihilation or a 
'higher twist effect' in lepton-hadron scattering, is expected to have the form of a 
series of inverse powers of the physical scale (see section 5). 

In practice most QCD calculations of observables are performed using finite- 
order perturbation theory, and calculations beyond leading order depend on the 
renormalisation scheme employed, implying a scheme-dependent strong-interaction 
scale A. It is conventional to work in the modified minimal subtraction scheme 



(MS scheme) [p^ , and to use the strong interaction scale Aj^s f° r fi ye active 
quark flavours. If one knows Aj^ one may calculate the strong coupling a s (Q 2 ) 
from the solution of the QCD renormalisation group equation [ |33|| : 

a ((J*) - 4n f 1 - ^ ln ( ln ^ 2 / A k)) , x (28) 
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Because of the large data samples taken in e + e~ annihilation at the Z° resonance, 
it has become conventional to use as a yardstick at s (M§), where Mz is the mass 
of the Z° boson; Mz ~ 91.2 GeV | 34j . Tests of QCD can therefore be quantified 



in terms of the consistency of the values of a s (M§) measured in different exper- 
iments. The 'QCD-challenged' reader may like to think of a s (M|) as being 'the 
sin 2 9w of strong interactions'. 

In e + e~ annihilation a s (M|) has been measured from inclusive observables 
relating to the Z° lineshape and to hadronic decays of the r lepton, as well 
as from jet-related hadronic event shape observables, and scaling violations in 
inclusive hadron fragmentation functions. 

4.2 R and the Z° Lineshape 

For the inclusive ratio R = <r(e + e~ — > hadrons) /cr(e + e~ — > fi + fi~), the SM 
electroweak contributions are well understood theoretically and the perturba- 
tive QCD series has been calculated up to O(a^) |35| for massless quarks and 



up to 0(a 2 ) including quark mass effects [P^] ; the large size of the O(af) term 
is potentially a cause for concern about the degree of convergence of the series. 
Closely-related observables at the Z° resonance are: 
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• the Z G total width, Yz 

• the pole cross section, cx° = \2ixT ee Thad/ M Z T 2 Z 

• the ratio of hadronic to leptonic Z° decay branching ratios Ri = Yhad/^u 
which all depend on the Z° hadronic width: 



T had = 1.671 (l + ai (^) + a 2 (^) 2 + a 3 (^) 3 + . . . J (29) 

where: a\ = 1, a 2 = 0.75 and = —15.3. In these cases the non-perturbative 
contributions are expected to be 0(1/ Mz) and are usually ignored. A concern 
is that recent measurements of observables that probe the electroweak couplings 
of the Z° to b and c quarks deviate slightly from SM expectations []3~7|| . Since 
these couplings must be known in order to extract a s (M z ), this effect, whatever 
its origin, is a potential source of bias P3"] . Further analysis is in progress from 
the SLC and LEP experiments and the situation is not yet resolved. 

Proceeding nonetheless, the procedure adopted [0 is to perform a global SM 



fit to a panoply of electroweak data that includes the W and top quark masses 
as well as the Z° observables relating to the lineshape, left-right production 
asymmetry, decay fermion forward-backward asymmetries, branching ratios to 
heavy quarks, and r polarisation. The free parameters are the Higgs mass, Mniggs, 
which contributes to X EW , and a s (M z ). Data presented at the 1996 summer 
conferences yield the results shown in Fig. 25 ||37|| , from which the positively- 
correlated results Mffiggs = 149^82° GeV and 

a s {M 2 z ) = 0.1202 ±0.0033 (exp.) (30) 

are obtained. The a s (M|) value is lower than the corresponding results presented 
at the 1995 conferences a s {M 2 z ) = 0.123±0.005, and at the 1994 conferences, 
o; s (M|) 0.125 ±0.005 |39|, whose large central values were partly responsible for a 



supposed discrepancy between 'low-Q 2 ' and 'high-Q 2 ' a s (M z ) measurements [|40 



The change between 1995 and 1996 is due to a combination of shifts in the values 
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of the Z° lineshape parameters, redetermined in light of the recalibration of the 
LEP beam energy due to the 'TGV effect' [37H, and a change in the central value 
of Mmggs at which a s (M§) is quoted, from 300 GeV (1995) to the fitted value 



149 GeV (1996). A detailed study of theoretical uncertainties implies [|T| that 
they contribute at a level substantially below ±0.001. Since data-taking at the 
Z° resonance has now been completed at the LEP collider the precision of this 
result is not expected to improve further. 




Figure 25: Results of a global fit of the Standard Model to electroweak observ- 
ables [jSTfl ; the 1- and 2-standard deviation contours are shown in the a s (M§) vs. 
M H iggs plane. 
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4.3 Hadronic r Decays 



An inclusive quantity similar to R is the ratio R T of hadronic to leptonic decay 
branching ratios, Bh and B\ respectively, of the r lepton: 

„ Bh 1 — Bp — B,, , . 

Rr = = 1 "~ (31) 

where B e and B^ can either be measured directly, or deduced from a measurement 
of the t lifetime t t . In addition, a family of observables known as 'spectral mo- 
ments' of the invariant mass-squared s of the hadronic system has been proposed 



where M r is the r mass. In this case the integrand can be measured independently 
of B e . It is easily seen that R T = 

i? r and R® have been calculated perturbatively up to O(af). However, because 
M T ~ 1 GeV one expects (eq. (28)) a s (M T ) ~ 0.3 and it is not a priori obvious 
that the perturbative calculation can be expected to be reliable, or that the non- 
perturbative contributions of 0(1/M T ) will be small. In recent years a large 
theoretical effort has been devoted to this subject; see eg. E3, [43], E4|. 



The ALEPH Collaboration derived R T from its measurements of B e , B^, and 
r T , and also measured the (10), (11), (12), and (13) spectral moments. A com- 
bined fit yielded |E§ a s (Mf ) = 0.124±0.0022±0.001, where the first error receives 



equal contributions from experiment and theory, and the second derives from un- 
certainties in evolving a s across the c and b thresholds. The OPAL Collaboration 
measured R T from B e , B^ and r T , and derived H] a s (M|) =0.1229^o;ool7 (exp.) 



iaorai (theor.). The CLEO Collaboration measured the same four spectral mo- 
ments as ALEPH and also derived R T using 1994 Particle Data Group values 
for B e , B^ and r T . A combined fit yielded [|7| a s (M|) = 0.114 ± 0.003. This 
central value is slightly lower than the ALEPH and OPAL values. If more re- 
cent world average values of B e and B^ are used CLEO obtains a higher central 

39 



a s (Mf ) value [g7|. Averaging the second CLEO result and the ALEPH and OPAL 
results by weighting with the experimental errors, assuming they are uncorrelated, 
yields: 

a s (M|) = 0.122 ±0.001 (exp.) ±0.002 (theor.). (33) 

This is nominally a very precise measurement, although recent studies have ruggested 
that additional theoretical uncertainties may be as large as ±0.006 



4.4 Hadronic Event Shape Observables 

As discussed in Section 3.2, in e + e~ annihilation the rate of 3-jet production: 

r 3 = (34) 

Chad 

is directly proportional to a s and can hence be used to determine a s . In order to 
make a meaningful measurement that can be compared with those just discussed 
one must calculate -R3 to at least next-to-leading order in a s , i.e. to O(a^). The 
relevant contributing Feynman diagrams are shown in Fig. 26; these form the 
basis of the O(a^) calculation of R3 BO, BO, . 



4.4.1 Definition of Jets and Event Shape Measures 

The task is, in principle, straightforward. One must count the number of 3-jet 
events and divide by the total number of hadronic events to obtain R3, then 
compare with the theoretical prediction to obtain a s . However, it is immediately 
apparent that one cannot simply define the jet multiplicity of events on the basis 
of a visual inspection! On the experimental side, the classic 'Mercedes-Benz' 3-jet 
event measured in a detector is rather rare; many events contain broad particle 
flows that might be classified as a single jet by one observer but as two or more jets 
by another observer. Moreover, in QCD the Bremsstrahlung spectrum of parton 
radiation peaks at small angles and is continuous. Hence even theoretically the 
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(a) 




Figure 26: Tree-level and one-loop Feynman diagrams that contribute to 3-jet 
observables up to 0(a 2 ) in QCD perturbation theory. 

issue of when a radiated parton is sufficiently energetic, and at a sufficiently wide 
angle relative to its parent, so as to be resolved as a separate jet is not without 
ambiguity. After due Cartesian deliberation one pragmatically concludes that one 
needs an algorithmic definition of a jet that can be applied to hadrons recorded in 
a detector, as well as to partons in perturbative QCD calculations, and a sensible 
recipe to translate between the two. 

A convenient solution is provided by iterative clustering algorithms in which 
a measure z/ij, such as invariant mass-squared/Q 2 , is calculated for all pairs of 
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particles i and j in an event, and the pair with the smallest yij is combined into 
a single 'particle'. This process is repeated until all pairs have exceeding a 
value y c , and the jet multiplicity of the event is defined as the number of particles 
remaining. For a sample of events the n-jet rate R n is then defined as the number of 
n-jet events divided by the total number of events. This number is not a constant, 
but rather depends on the choice of algorithm and on the y c value. The y c - 



dependence is illustrated in Fig. 27 for jets defined using the JADE algorithm [52 



applied to SLD data One can think of y c as the 'jet resolution' scale. Large y c 



values correspond to poor eyesight, most events look 2-jet-like, and hence R 2 ~ 1. 
Small y c values correspond to good eyesight, a richer jet structure is discernible, 
and -R3 and -R4 are non-zero. It should be noted, however, that from an operational 
point-of-view the data points shown in Fig. 27 are awkward to handle in that they 
are correlated between different y c values. A more convenient observable is the 
differential 2-jet rate: 

D 2 (y c ) = (MVc) ~ MVc - Ay c ))/Ay c (35) 

which is a measure of the rate of events that change their classification between 
2-jet-like and > 3-jet-like as y c is varied across the range Ay c . D 2 (y c ) is illustrated 
in Fig. 28. 



In fact several variations of the JADE algorithm have been suggested |p4 
these differ in the definition of the resolution measure y^, and/or in the 'recom- 
bination scheme' prescription for combining two particles that are unresolvable. 
A full discussion is beyond the scope of these lectures, but it is important to note 
that the 'E', 'EO', 'P' and 'PO' variations of the JADE algorithm, as well as the 
'Durham' ('D') and 'Geneva' ('G') algorithms, are all collinear- and infra-red- 
safe observables, which, for our purposes, means that they can be calculated in 
perturbative QCD |p3 



More generally one can define other infra-red- and collinear-safe measures of 
the topology of hadronic final states; a list of 15 such observables is given in 
Table 3. Thrust has already been encountered in Section 3 and is related to the 
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Figure 27: Dependence on the jet resolution parameter y c of the n-jet rates R Ti 



measured using SLD data [53] with the JADE algorithm 



longitudinal momentum flow in events: 

1 = max — — , (3o) 





p 


i-n T \ 




Pi 





where pi is the momentum vector of particle i, and n-r is the thrust axis to be 
determined. It is useful to define r = 1 — T. For back-to-back two-parton final 
states r is zero, while < r < | for planar three-parton final states. Spherical 
events have r = |. An be found to maximize the momentum sum 

transverse to n^, and defined to be perpendicular to the two axes 

nj- and n ma j. The variables thrust-major T ma j and thrust-minor T min are obtained 
by replacing Ht in Eq. (|36| ) by n ma j or n m j n , respectively. The oblateness O is 
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Figure 28: Dependence on the jet resolution parameter y c of the differential 2-jet 
rate D 2 [E| . 



then defined by [gC 



O T ma j T m i n . (37) 
Other measures are related to jet masses, and energy-energy correlations between 



particles; for a discussion see eg. [[57 . 



The observables are all constructed to be directly proportional to a s at lead- 
ing order, and so are potentially sensitive measures of the strong coupling. The 
0(«g) QCD prediction for each of these observables X can be written |5T| : 

1 d(J ( TD/V\f a s\ 2 



a dX 
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Observable 


symbol 


1 - Thrust 


T 


Heavy jet mass 


P 


Jet broadening: 




Total 


Bj> 


Wide 


B w 


Oblateness 





C-parameter 


C 


Differential jet rates: 






E 




EO 




P 




PO 




D 




G 


Energy-energy correlations 


EEC 


Asymmetry of EEC 


AEEC 


Jet cone energy fraction 


JCEF 



Table 3: Fifteen infra-red- and collinear-safe measures of the topology of 
e + e~ hadronic final states. 

so that a s can be determined from each. Though these observables are intrinsically 
highly correlated, by using all 15 to study a s one is attempting to maximise the 
use of the information in complicated multi-hadron events, and in some sense is 
making a more demanding test of QCD than by using only one or two observables. 
Moreover, it will be seen that the study of many observables is essential, as it 
may expose systematic effects. Finally, the a s determination from hadronic event 
shape observables is based on the information content within 3-jet-like events, and 
is essentially uncorrelated with the measurements from the Z° lineshape which 
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are based on event- counting of predominantly 2-jet-like final states. 

The technology of this approach has been developed over the past 15 years 
of analysis at the PETRA, PEP, TRISTAN, SLC and LEP colliders, so that the 
method is considered to be well understood both experimentally and theoretically. 
Note, however, that before they can be compared with perturbative QCD predic- 
tions, it is necessary to correct the measured distributions for any bias effects 
originating from the detector acceptance, resolution, and inefficiency, as well as 
for the effects of initial-state radiation and hadronisation, to yield 'parton-level' 
distributions. 




isation 
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Figure 29: Schematic of hadron production in e + e annihilation. 
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4.4.2 Hadronisation and Monte Carlo Models 



A schematic of hadron production in e + e annihilation is shown in Fig. 29. One 
may divide this process into several phases: 

1. A hard electroweak process in which the primary quark and antiquark may be 
produced off mass-shell: 

e + e~ — ► qq 

2. Perturbative QCD evolution of the primary qq via parton Bremsstrahlung: 
qq — > several q, q, g 

3. Hadronisation of partonic system: 

(q, q, g)s — > primary resonances 

4. Decays of primary resonances into 'stable' particles: 

B, C, K°„ <f>, A, p . . . - 7r± K± p,p . . . (e± fi±, r± v) 

Phases 1 and 2 are generally agreed to be calculable 'respectably' using perturba- 
tive techniques applied to the electroweak theory and QCD, respectively. Phases 
3 and 4 are more problematic in that they are intrinsically non-perturbative pro- 
cesses that cannot in general be calculated from first principles. In the absence of 
non-perturbative calculations we are forced to rely on phenomenological models. 

Since it is also necessary in phase 4 to simulate the interaction of particles 
with detectors, which can only be done in a deterministic fashion, Monte Carlo 
event generators have been developed for the complete simulation of hadronic 
event production in e + e~ annihilation and are now essential components of data 
analysis. I shall discuss only the two most widely used generators JETSET |58 



and HERWIG |p9jl ; other generators are described in [pOj] , and will be discussed 
later by Buchanan [5T|. I shall not discuss at all the GEANT program [JS2[| , which 



is widely used for the simulation of the geometry and material response of particle 
detectors. The philosophy here is to outline the main features of these generators 
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in the context of their use as tools in understanding and correcting the data; no 
attempt will be made to justify these models on phenomenological grounds, and 
the outline will necessarily be brief. 

Both JETSET and HERWIG implement electroweak matrix elements for the 
production of a primary qq, as well as a perturbative QCD 'parton shower' evolu- 
tion of the system into a set of low- virtual-mass quarks and gluons. More formally, 
the latter is based on a probabilistic parton branching process that is derived from 
a leading + partial next-to-leading logarithmic resummation of the QCD matrix 
elements []63] . JETSET and HERWIG implement the parton branching process 
slightly differently, a discussion of which is beyond the scope of these lectures, but 
both generators have a parameter A that characterises the scale of strong inter- 
actions, as well as a parameter Q that characterises the minimum virtual-mass 
scale of the parton evolution. 

A schematic of the hadronisation process as implemented in HERWIG is shown 
in Fig. 30. At the termination of the parton shower pairs of partons are associateed 
into colourless clusters; these then undergo phase-space decay to produce stable 
pions, kaons and baryons. Clusters with mass larger than a parameter M c ; are 
split into two before the phase-space decay. Additional parameters control the 
properties of heavy (B or C) hadron decay |59|] . 

JETSET implements the 'Lund string model' of jet fragmentation pi], illus- 
trated in Fig. 31. In this case the colour field between partons at the end of 
the parton shower is represented as a one-dimensional massless relativistic string. 
String pieces terminate at quarks and antiquarks, and gluons are represented by 
momentum-carrying 'kinks' in the string. The string is fragmented iteratively 
according to the recipe: 

/(*) oc l(l-*)»exp(-6^) (39) 

where z is the fraction of the quantity E + pii of a parent string piece taken by 



the daughter, m±_ = Jp\ + m 2 , '_!_' and ' ||' refer to the string axis, and a and 
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Figure 30: Schematic of hadronisation in HERWIG. 

b are parameters. Momentum transverse to the string axis, p±, is introduced in 
an ad hoc fashion using a Gaussian probability distribution. A large number of 
additional parameters is used to fine-tune the relative production of particles such 
as strange, pseudoscalar, and vector mesons, as well as strange and non-strange, 
octet and decuplet baryons |58 . 



4.4.3 Data Correction 

For the a s analysis we wish to use these event generators to understand the effect 
of the hadronisation process on the hadronic momentum flow in events, and to 
correct for any bias, as well as to understand the influence of the response of 
the detector. One conventional approach involves using a sample of Monte Carlo 
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Figure 31: Schematic of hadronisation in JETSET. 

events to calculate bin-by-bin correction factors, and then applying these to the 
measured distribution. For a distribution D(X), the correction for detector effects 
is defined: 

where HAD and DET refer to the simulated distribution at the hadron-level 
and detector-level phases, respectively. The correction for hadronisation effects is 
analogously defined: 

nMC f V \ _ Dpart{.X) , . 

UhadK^) 

where PART refers to the simulated distribution at the parton-level. The data 
distribution corrected back to the parton-level is given by: 

D Data \X) = C^ D (X)-C^ T (X)-D Data (X), (42) 

and D Data (X) can be compared with perturbative QCD. More sophisticated cor- 
rection procedures can also be defined; see eg. |p5| . 
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As with any correction procedure one must take care not to introduce bias 
from implicit model-dependence, and must estimate the systematic uncertainties 
involved. A prerequisite is that the simulation describe the distribution measured 
in the detector! The parameters of the detector simulation, as well as of the event 
generator itself, should then be varied, the stability of the correction factors exam- 
ined, and systematic errors assigned accordingly. An example of a raw measured 
D2 distribution from SLD ||57|| , compared with simulations based on JETSET and 
HERWIG, is shown in Fig. 32. The corresponding corrected distribution, and the 
correction factors, are shown in Fig. 33. 



10 1 

D 

a 

^10° 
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0.1 0.2 0.3 
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Figure 32: Measured D2 distribution |57| compared with JETSET and HERWIG 
predictions. 
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Figure 33: (a) Corrected D 2 distribution ||57l ; (b) hadronisation correction factor; 
(c) detector effects correction factor. In (a) the line shows a fit to O(a^) QCD. 

4.4.4 Comparison with Perturbative QCD 

A fit of O(ag) perturbative QCD to the D 2 distribution is shown in Fig. 33; it 
yields a s {M\) = 0.1175 ±0.0007 (stat.) ±0.0027 (syst.) CT. One can repeat 



this procedure for all 15 observables listed in Table 3 and derive in each case a 
fitted value of a s (M§); these are shown in Fig. 34a. The distressing result of 
this exercise is that the a s (M§) values so determined are not internally consistent 
with one another! A measure of the scatter among the results is given by the 
r.m.s. deviation of ±0.008, which is much larger than the experimental error of 
±0.003 on a typical observable. An exciting, though remote, possibility is that 
we have observed a spectacular breakdown of QCD! A more likely explanation is 
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that some systematic effect that we have not yet considered is at work. In fact 
an implicit assumption was made in deriving the results shown in Fig. 34a that 
relates to the arcane issue of choosing the renormalisation scale in QCD. 
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Figure 34: (a) Values of a s (M|) determined |57j by fitting O(a^) QCD predic- 
tions to 15 hadronic event shape observables using a fixed value of the renormal- 
isation scale ji — Q: the results are clearly inconsistent within the experimental 
errors, (b) Renormalisation scale uncertainties. 



4.4.5 Renormalisation Scale Uncertainty 

For any observable, truncation of the QCD perturbation series at finite order 
causes a residual dependence on the (scheme-dependent) renormalisation scale 
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/i. This parameter is formally unphysical and should not enter at all into an 
exact infinite-order calculation, and its value is arbitrary For the event shape 
observables an explicit /i-dependence enters the next-to- leading coefficient: 



1 da 



A(X) 



B{X) + A(X)2irb V/Q 5 



(43) 



a dX " s " y \ 2tt J ' v """ / ■•-»—'-/-« >y 2 tt 

so that a measurement of a s must be in the context of some chosen value of \i. 
This is illustrated in Fig. 35, where the value of from fits to D 2 is shown 
as a function of the choice of /i; there is clearly a strong /z-dependence. The 
top portion of Fig. 35 shows the corresponding Xdof f° r each fit; amusingly the 
data show no preference for any particular value of \i provided it is larger than 
V0.001Q 2 - A full discussion of the form of the /i-dependence is beyond the scope 



of these lectures; see [p6| . Figures of the /i-dependence for the other observables 
can be found in I57 . 
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Figure 35: Dependence of (a) Xdof an d (b) ^-Ws 011 ^ n e value of \i chosen in fits 
to the SLD D 2 distribution 
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A consensus has arisen among experimentalists that the effect of missing 
higher-order terms can hence be estimated from the dependence of a s (M§) on 
the value of // assumed in fits of the calculations to the data, and a renormali- 
sation scale uncertainty is often quoted. This procedure, well-motivated in that 
the /x-dependence caused by the truncation of the perturbation series would be 
cancelled by addition of the higher-order terms, is, however, arbitrary, and is not 
equivalent to knowledge of the size of the a priori unknown terms. In cases where 
scale uncertainties are considered this arbitrariness is manifested in the wide vari- 
ation among the ranges and central values of fi chosen by different experimental 



groups, see eg. []67[]; in other cases this source of uncertainty is not included in 
the errors. Different a s (M|) results with similar experimental precision can hence 
be quoted with different total errors depending on the procedure adopted for as- 
signing the theoretical uncertainties. The interpretation of the central values and 
errors on a s (M|) measurements is hence not always straightforward. The SLD 
estimate of the renormalisation scale uncertainty for each observable is shown in 
Fig. 34b. It is apparent that the scale uncertainty is much larger than the experi- 
mental error, and that the a s (M|) values are consistent within these uncertainties. 
Though this is comforting, in that it indicates that QCD is self-consistent, the nec- 
essary addition of large theoretical uncertainties to otherwise precise experimental 
measurements is frustrating, at least to experimentalists! 

The best resolution of the scale ambiguity would be to reduce its effect by 
calculating observables to higher order in perturbation theory. Though this is 
in principle possible, the large number of Feynman diagrams involved renders 
the task difficult and unattractive. In e + e~ annihilation only the i?-related ob- 
servables and the r hadronic decay ratio R T , have been calculated exactly up to 
0(al). For the hadronic event shape observables 0(al) contributions have not yet 
been calculated completely. However, for six observables (indicated in Table 3) 
improved calculations can be formulated that incorporate the resummation |)8| 
of leading and next-to- leading logarithmic terms matched to the 0(ai?) results. 
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The matched calculations are expected a priori both to describe the data in a 
larger region of phase space than the fixed-order results, and to yield a reduced 
dependence of a s on the renormalization scale. This is illustrated in Fig. 36 for 
the case of thrust (r). Though not well described by the 0(a 2 s ) calculation, the 
low-r region is well reproduced when resummed contributions are included. 

Application of other approaches to circumvent the scale ambiguity in a s mea- 
surement, involving the use of 'optimised' perturbation theory' and Pade 
Approximants |7(J , can be found in |J7|, [7]J respectively. 



i i i i i i r 

SLD C>(a!) 

Resummed + O(o? s ) 

mod. InR 
• Data 




Figure 36: Illustration of the need for resummed contributions: the r distribution 
measured by SLD ]57|| . At low r the O(a^) calculation is unable to describe the 
data unless resummed terms are considered. 
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4.4.6 Summary of a s Measurements 



Hinchliffe has reviewed the various hadronic event shapes-based measurements 
from experiments performed in the cm. energy range 10 < Q < 91 GeV, utilising 
both O(aig) and resummed calculations, and quotes an average value of a s (M|) = 



0.122 ±0.007 [33], where the large error is dominated by the renormalisation scale 
uncertainty, which far exceeds the experimental error of about ±0.002. Schmelling 
has also compiled the measurements, including the recent results from the LEP-II 
run at Q ~ 133 GeV |72]] , and quotes a global average []73[ a s {M\ ) = 0.121±0.005, 
in agreement with |j33"l, but assuming a more aggressive scale uncertainty. 



4.5 Scaling Violations in Fragmentation Functions 

Though distributions of final-state hadrons are not, in general, calculable in per- 
turbative QCD, the Q 2 -evolution of the scaled energy (x p = 2E/Q) distributions 
of hadrons, or 'fragmentation functions', can be calculated and used to determine 
a s . In addition to the usual renormalisation scale /j,, a factorisation scale fiF must 
be defined that delineates the boundary between the calculable perturbative, and 
incalculable non-perturbative, domains. Additional complications arise from the 
changing composition of the underlying event flavour with Q due to the differ- 
ent Q-dependence of the 7 and Z° exchange processes. Since B and C hadrons 
typically carry a large fraction of the beam momentum, and contribute a large 
multiplicity from their decays, it is necessary to consider the scaling violations 
separately in b, c, and light quark events, as well as in gluon jet fragmentation. 

In an early analysis |73j] the DELPHI Collaboration parametrised the fragmen- 
tation functions using the O(ct^) matrix elements and the string fragmentation 
model implemented in JETSET [Q. They fitted data in the range 14 < Q < 91 
GeV to determine a s (M|) = 0.118 ±0.005, where the error is dominated by vary- 
ing \x in the range 0.1 < ji/Q < 1. The ALEPH Collaboration used its Z° data to 
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constrain flavour- dependent effects by tagging event samples enriched in light, c, 
and b quarks, as well as a sample of gluon jets [[75]. The fragmentation functions 
for the different flavours and the gluon were parametrised at a reference energy, 
evolved with Q according to the perturbative DGLAP formalism calculated at 
next-to-leading order [f?6|], in conjunction with a parametrisation proportional 
to 1/Q to represent non-perturbative effects (Section 5), and fitted to data in 
the range 22 < Q < 91 GeV (Fig. 37). They derived a s (Mf ) = 0.126 ± 0.007 
(exp.) ±0.006 (theor.), where the theoretical uncertainty is dominated by varia- 
tion of the factorisation scale /ip in the range — 1 < lnfi 2 F /Q 2 < 1; variation of 
the renormalisation scale in the same range contributed only ±0.002. DELPHI 
has recently reported a similar analysis [77] yielding a s {M\) = 0.1211;q;qq7 (exp.) 
±0.010 (theor.). Curiously, although a similar range as ALEPH, 0.3 < [i/Q < 3, 
was used to examine variation of the renormalisation and factorisation scales, here 
the renormalisation scale dominates the theoretical uncertainty, with a contribu- 
tion of ±0.009, in contrast to ±0.002 from factorisation. Combining the ALEPH 



and later DELPHI results, assuming uncorrelated experimental errors, yields [ plf : 
a s (M 2 z ) = 0.124 ± 0.005(exp.) ±0.010(theor.) (44) 



4.6 Comparison with Other Measurements of a s (M§) 

A summary of world a s measurements, all evolved to Q — Mz, is shown in 
Fig. |38] [|3T|. These are drawn from lepton-hadron scattering, hadron-hadron col- 
lisions, heavy quarkonia decays and lattice gauge theory, as well as e + e~ annihila- 
tion. In addition to being relatively precise, the e + e~ results have the invaluable 
feature that they bracket the Q-range of the experiments, from around 1 GeV for 
r decays to around 100 GeV for Z° production, providing the largest lever-arm 
for tests of consistency of a s (M|) measured at different energy scales. It is clear 
that, within the uncertainties, all results are consistent with one another. 
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Figure 38: Summary of world a s (M|) measurements | 3l| . The results are ordered 
vertically in terms of the hard scale Q of the experiment. 
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Taking an average over all 17 measurements assuming they are independent, by 
weighting each by its total error, yields a s (M|) = 0.118 with a x 2 of 6.4; the low x 2 
value reflects the fact that most of the measurements are theoretical-systematics- 
limited. Taking an unweighted average, which in some sense corresponds to the 
assumption that all 17 measurements are completely correlated, yields the same 
result. The r.m.s. deviation of the 17 measurements w.r.t. the average value 
characterises the dispersion, and is ±0.005. In a quantitative sense, therefore, 
QCD has been tested to a level of about 5%. 

If further progress is to be made in testing QCD, future measurements of 
a s (M|) should aim for substantially improved precision. The prospects for achiev- 
ing 1%-level measurements are discussed in detail elsewhere [J78| ]. Lattice QCD 
determinations may reach this precision within the next few years. A precise 
a s (M§) measurement has yet to emerge from the TeVatron, but feasibility stud- 
ies are in progress and appear promising. Deep-inelastic scattering and e + e~ 
annihilation will probably require higher-energy facilities, as well as significant 
theoretical effort with regard to O(a^) perturbative contributions. An a s (M|) 
measurement at a high-energy e + e~ collider will be discussed in Section 6.2. 



5. Towards a Theory of Hadronisation 

We expect that the strong coupling becomes large in long-distance (low-Q) q— q 
interactions such that finite-order perturbation theory is no longer valid. Lattice 
gauge theory |79[ is the only practical non-perturbative calculational tool available 
today. It is presently limited in applicability to static properties of hadrons, such 
as masses and decay constants, although in principle it might eventually be applied 
to the dynamical process of hadronisation. 

From the operator product expansion (OPE) one expects (see eg. HSOfl ) that 
the expectation value of an observable O may be written: 

<0> = + E | (45) 
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During the past 15 years much theoretical effort has been focussed on the pertur- 
bative component represented by the first term in this equation. More recently 
attention has turned to the 'power corrections' represented by the second term, 
whose origin is intrinsically non-perturbative. In particular, attempts have been 
made to evaluate power corrections for e + e~ observables. An illustration of the 
potential of such an approach is provided by Fig. |39] The ad hoc addition 
of a 1/Q term to the 0(a;?) QCD prediction describes the energy-dependence 
of < 1 — T > remarkably well. It will be seen in Section 6.2 that the inverse 
power-law behaviour of hadronisation effects has important consequences for a 
precise a s measurement at a high energy e + e~ collider. The explicit calculation 
of leading power corrections [BD| hence represents our first tentative step towards 
a consistent theoretical treatment of hadronisation. 



6. QCD at a High Energy e + e Collider 
6.1 Introduction 

Since QCD is our theory of strong interactions it would be irresponsible not to 
test it at the highest energy scales available in different hard scattering processes. 
For this reason testing QCD at a 0.5-1.5 TeV e + e~ collider ('XLC') is mandatory. 



For a detailed discussion see 82 



Precise determination of the strong coupling a s is key to a better understand- 
ing of high energy physics. The current precision of a s (M|) measurements, limited 
to about 5% (Section 4.6), results in the dominant uncertainty on our prediction 
of the energy scale at which grand unification of the strong, weak and electro- 
magnetic forces takes place. An a s (M|) measurement of 1% precision may be 
possible at a high energy e + e~ collider. Such a measurement would also allow 
improved determination of the mass and width of the top quark from the thresh- 
old behaviour of the tt cross-section. Measurements of hadronic event properties 
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Figure 39: Energy- dependence of the mean value of 1 — T. The ad hoc addition 
of a 1/Q term to the O(a^) QCD prediction describes the data [|3T |. 



at high energies, combined with existing lower energy data, would allow one to 
test further the gauge structure of QCD by searching for anomalous 'running' of 
observables, such as the rate of production of events containing three jets, and to 
set limits on models which predict such effects, for example those involving light 
gluinos which are difficult to exclude by other means. 

Gluon radiation in tt events is expected to be strongly regulated by the large 
mass and width of the top quark; ttg events will hence provide an exciting new 
domain for QCD studies. As a corollary, measurements of gluon radiation patterns 
in ttg events may provide valuable additional constraints on the top quark decay 
width. Furthermore, searches could be made for anomalous chromo-electric and 
chromo- magnetic moments of quarks ||83|| , which effectively modify the rate and 



63 



pattern of gluon radiation, and for which the phase space increases as the cm. 
energy is raised. Finally, polarised electron beams will be exploited at high energy 
e + e~ colliders and will allow tests of symmetries using multi-jet final states |84 . 



6.2 Is a 1%— level Measurement of a s (M|) Possible? 

It is interesting to consider whether a measurement of a s (M|) at the 1%-level 
of precision is possible at the XLC Consider the SLD a s (M|) measurement, 
discussed in Section 4.4, based on 15 hadronic event shape observables measured 



with a data sample comprising approximately 50,000 hadronic events [57| : 



a s {M 2 z ) = 0.1200 ± 0.0025 (exp.) ± 0.0078 (theor.) (46) 

where the experimental error is composed of statistical and systematic compo- 
nents of about ±0.001 and ±0.002 respectively, and the theoretical uncertainty 
has components of ±0.003 and ±0.007 arising from hadronisation and missing 
higher order terms, respectively. Now consider 'scaling' this result to estimate the 
precision of a similar measurement at Q = 500 GeV. 

• Statistical error: At design luminosity the 500 GeV XLC would deliver roughly 
100,000 qq (q=u,d,s,c,b) events per year (Section 6.4), implying that a statistical 
error on a, (Mf) well below ± 0.001 could be obtained. 

• Systematic error: This results primarily from the uncertainty in modelling 
the jet resolution of the detector. The situation may be improved at the XLC by a 
combination of building better detectors and benefitting from improved calorime- 
ter energy resolution for higher energy jets. It is not unreasonable to suppose that 
the current systematic error of roughly ±0.002 could be reduced by a factor of 
two. 

• Hadronisation uncertainty: From the discussion in Section 5 it can be seen 
that non-perturbative corrections to jet final states in e + e~ annihilation can be 
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parametrised in terms of inverse powers of the hard scale Q. At leading order, 
perturbative evolution is proportional to 1/lnQ. Hence for a generic observable 
X the ratio of non-perturbative to perturbative QCD contributions is dominated 
by a term of the form: 

5X non-pert j n Q 

~ — — . (47) 

Xpert Q V ) 

Increasing Q from 91 GeV to 500 GeV causes this ratio to decrease by a factor of 5, 
implying that hadronisation corrections in the '3-jet region' of observables should 
be of order 2% at XLC The conclusion of this analysis is reinforced by explicit 



simulation of hadronisation effects, illustrated in Fig. |40] |85j for thrust. Assuming 
that these corrections can be estimated to better than ±50%, the hadronisation 
uncertainty should contribute less than 1% to the error on a s (M§). 

• Uncertainty due to missing higher orders: Currently perturbative QCD 
calculations of hadronic event shapes are available complete up to 0(a 2 s ). Since 
the data contain knowledge of all orders one must estimate the possible bias inher- 
ent in measuring a s (M|) using the truncated QCD series (Section 4.4.5). Since 
the missing perturbative terms are O(a^), and since at Q = 500 GeV a s is ex- 
pected to be about 25% smaller than its value at the Z°, one naively expects the 
uncalculated terms to be almost a factor of two smaller at the higher energy, lead- 
ing to an estimated uncertainty of ±0.004 on a s (500 GeV). However, translating 
to the yardstick a s (M|) yields an uncertainty of ±0.006, only slightly reduced 
compared with the current uncertainty. 

From this simple analysis it seems reasonable to conclude that achievement 
of the luminosity necessary for 'discovery potential' at the XLC will result in a 
qq event sample of sufficient size to measure a s (M|) with a statistical uncertainty 
of better than 1%. Construction of detectors superior in performance to those in 
operation today at SLC and LEP may be necessary in order to reduce systematic 
errors to the 1% level. Hadronisation effects should be significantly smaller, im- 
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Figure 40: Estimate of the hadronisation correction factor (Section 4.4.2) for 



thrust at Q = 91 GeV and 500 GeV At 500 GeV the factor barely deviates 
from unity for most of the kinematic range. 

plying a sub-1% uncertainty However, unless O(af) contributions are calculated, 
Q! S (M|) measurements at 500 GeV will be limited by theoretical uncertainties to 
a precision of ±0.006, only marginally better than that achieved at present. 

6.3 Top Quark Mass Determination and a s 



It is clear that the value of a s controls the shape of the strong potential that binds 
quarkonia resonances. In the case of tt production near threshold, the large top 
mass m t , and hence large decay width T, ensure that the top quarks decay in a 
time comparable with the classical period of rotation of the bound system, making 
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the toponium resonance a very short-lived phenomenon, and washing out most of 
the resonant structure in the cross-section. The shape of the tt cross-section near 
threshold hence depends strongly not only on the top mass, but also on a s . 
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Figure 41: Simulation of a simultaneous measurement of a s (M|) and m t at a 
high-energy e + e~ collider |p6[ : (a) tt production cross section; (b) correlation 



between fitted values. 



Fits to simulations of measurements of this cross-section have shown |86[ that 
the top mass so determined is strongly correlated with the assumed value of 
a s {M\\ This is illustrated in Fig. |41]. The European Top Quark Working Group 
has updated these simulations for the latest measured values of the top mass and 
has shown |87j that a simultaneous determination of m t and a s (M|) by fitting 
to the threshold cross-section measured with one design- year of luminosity yields 
statistical precisions of ±250 MeV/c 2 and ±0.006 on m t and a s (M|), respectively 
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Fixing a s (M§) to 0.120 reduces the error on m t by a factor of 2. Since this tech- 
nique would yield a measurement of a s (M§) no more precise than those made 
today, and since systematic uncertainties may be large and have not yet been 
considered, a more sensible strategy would be to measure a s (M|) as precisely as 
possible, as described in the previous section, and to use this value to allow better 
determination of the top quark parameters. 



6.4 Energy Evolution Studies 

The non-Abelian gauge structure of QCD implies that as the hard scattering scale 
Q increases, the strong coupling decreases roughly as 1/lnQ. Existing hadronic 
final states data from e+e" annihilation at the PETRA, PEP, TRISTAN, SLC 
and LEP colliders span the range 14 < Q < 170 GeV, although hadronisation 
uncertainties are large on the data below 25 GeV. A 1.5 TeV e + e~ collider would 
increase the lever-arm in 1/hiQ by almost a factor of two, hence allowing detailed 
study of the energy evolution of QCD observables that are proportional to a s , 
such as the rate of production of final states containing three hadronic jets, R3. 
This would provide not only a test of the fundamental structure of SU(3) QCD, 
but also a search-ground for new physics that might produce 'anomalous' running. 

One such possibility is the existence of a light, electrically neutral, coloured 
fermion that couples to gluons, often called a 'light gluino' and denoted by g. 
The existence of such a particle would manifest itself via a modification of gluon 
vacuum polarisation contributions involving fermion loops, effectively increasing 
the number of light fermions entering into the QCD /3-function. At one-loop level 
the effective number of flavours would change from Np to Np + 3A^, where Ng 
is the number of families of light gluinos, causing a decrease in the running of 
a s as a function of Q. The existence of a light gluino of mass between 2 and 5 



GeV/c has not been excluded by searches with current data [85|. A simulated 



measurement of R 3 at Q = 500 GeV, corresponding to one design-luminosity- 
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year, is shown in Fig. [|2] |85|], together with existing measurements, plotted as a 
function of 1/lnQ. The presence of one family of light gluinos of mass 2 GeV/c 2 
would cause an increase in the predicted value of R3 at 500 GeV by 10%. A 1%- 
level measurement of R 3 , as discussed in the previous section, would allow this 
difference to be measured with a significance of many standard deviations. 
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Figure 42: Energy-evolution of the 3-jet rate R3 f85[. For illustration an 



0(«g) QCD fit, as well as a fit allowing the possibility of one family of light 
gluinos, is shown. The simulated data point at Q = 500 GeV would add consid- 
erable lever-arm. 

It should be noted, however, that data from a number of experiments at differ- 
ent e + e~ colliders contribute to Fig. f|2[ Some of these data were recorded more 
than 10 years ago, were treated differently by the various experimental groups, 
and have relatively large systematic errors that are at least partly uncorrelated 
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from point to point. Furthermore, the sophistication and performance of parti- 
cle detectors constructed in the last decade has improved significantly, and it is 
reasonable to assume that future detectors will be even better. In addition, our 
understanding of the modelling of hadronisation effects and theoretical uncertain- 
ties has improved enormously as a result of studies at the Z°. Therefore, the 
precision of searches for anomalous running of QCD observables at XLC would 
be improved significantly if new data were taken at the lower cm. energies with 
the same detector and analysis procedures. 

In fact, if the luminosity of the 500 GeV XLC could be preserved at lower cm. 
energies, very large data samples would be recorded. Table |] []85| shows the num- 
ber of qq events delivered per day at various cm. energies by the XLC operating 
at the design luminosity of 5 x 10 33 cm~ 2 s -1 . At each energy more luminosity 
would be delivered per day than was recorded in total by the original dedicated 
colliders! This argument is of course naive, in that a collider designed to operate 
at a luminosity of 5 x 10 33 cm _2 s _1 at 500 GeV would not automatically be oper- 
able at the same luminosity at energies a factor of 5 or 10 lower; such capability 
would have to be designed from the outset. Furthermore, the requirements on 
the triggering and data processing capabilities of the detector are extreme by the 
standards of e + e~ annihilation, and this would also have to be designed from the 
start. Nevertheless, the prospect of running the XLC at the Z° resonance, or at 
even lower energies, for QCD studies, not to mention high-statistics electroweak 
physics measurements, is very attractive. 

6.5 Gluon Radiation in tt Events 

The large mass and decay width of the top quark serve to make the study of 
gluon radiation in tt events a new arena for testing QCD. The large mass acts as 
a cutoff for collinear gluon radiation, and the large decay width acts as a cutoff 
for soft gluon radiation, allowing reliable perturbative QCD calculations to be 
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cm. energy Q (GeV) 


qq events/day 


500 


1750 


91 


20,000,000 


60 


75,000 


35 


150,000 



Table 4: Number of qq events per day delivered by an e + e collider operating at 
a luminosity of 5 x 10 33 cm" 2 s _1 . 

performed; these effects are of course correlated. The latter case is particularly 
interesting. If the top width were infinite, top quarks would decay immediately 
to bottom quarks, and any gluons would be radiated from the secondary b's. If 
the top width were zero, top quarks would live forever and all radiation would 
be from the primary t's. In the case of a large but finite width, expected to be 
around 2 GeV for a top mass of 180 GeV/c 2 , gluon radiation in tt events will be 
a coherent sum of contributions from these two limiting cases, with a degree of 
coherence regulated by the top width itself. 

A theoretical study of tt production above threshold, assuming m t = 175 
GeV/c 2 at Q = 1 TeV, is illustrated in Fig. |43] [B8fl. This shows the angular distri- 
bution of 5 GeV gluons w.r.t. the tt axis for the kinematic configuration in which 
the decay b-quark travels backwards w.r.t. the t flight direction. The dependence 
of the radiation pattern on the top decay width is strong. Similar effects are 



predicted in the spectrum of gluon radiation in tt events around threshold [p9 
Measurement of such effects would yield not only a dramatic demonstration of 
quantum interference in strong interactions, but might also provide an essential 
cross-check on the value of the top quark decay width, which may prove difficult 
to disentangle from measurements of the tt threshold cross-section and top mo- 
mentum distributions, which also depend on a s and m t (section 6.3), as well as 
on the beam energy distribution. 
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Figure 43: Angular distribution of 5 GeV gluons w.r.t. tt axis, at Q — 1 TeV, 



illustrating the dependence on the top width T [|8S 
7. Concluding Remarks 



We have seen that e + e~ annihilation is an ideal laboratory for precise studies of 
QCD. One observes jets indicating the primary production of quarks and gluons, 
and one can measure precisely the quark and gluon spins. Multijet events allow 
the very gauge structure of QCD to be tested via measurement of the Casimir 
factors N c , T F , and Cp, leading us to the conclusion that QCD is the theory 
of strong interactions. One can then measure the single parameter of QCD, the 
coupling a s , from inclusive observables such as R, or equivalently the Z° line- 
shape parameters, and from hadronic r decays, as well as from event shape mea- 
sures and scaling violations in inclusive single-particle fragmentation functions. 
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These a s (M|) measurements are internally consistent, and agree with results 
from lepton-nucleon scattering, hadron-hadron collisions, and lattice gauge the- 
ory determined across a wide range of energy scales. 

There was no time to cover many interesting topics, including: differences 
between quark and gluon jets, tests of the flavour-independence of strong interac- 
tions, polarisation phenomena, particle multiplicities and correlations, production 
of B and C mesons and baryons, and production of identified hadrons such as 7r ± , 
K°, K ± , p/p, A, cf), K* etc. Some of these topics are discussed in other contributions 
to these proceedings |Tl], 90 1. 



Looking towards the future, tests of QCD will provide an important component 
of the physics programme at a future high energy e + e~ collider operating in the 
cm. energy range 0.5 < Q < 1.5 TeV. Measurement of a s (M|) at the 1% level of 
precision appears feasible experimentally, but will require considerable theoretical 
effort to calculate O(a^) contributions in QCD perturbation theory. A search for 
anomalous running of a s (Q 2 ), by operating the collider at different cm. energies, 
is an attractive prospect. Quantum coherence is expected to give rise to interesting 
gluon radiation patterns in tt events, which could be used to constrain the top 
quark decay width, and measurement of the gluon radiation spectrum would also 
constrain anomalous top quark chromomagnetic couplings. 

More immediately, the next generation of low energy e + e~ colliders, known as 
B factories, also has the potential to make a precise a s measurement from the 
i?-ratio at Q ~ 10 GeV, as well as from hadronic r decays. Even more precise 
tests of QCD in e + e~ annihilation will hence continue to enhance our confidence 
in the theory, and may even yield surprises. . . 
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